NASA Contractor Report 3968 


i 


NASA-CR-3968 19860013110 


Extensions to PIFCGT: Multirate 
Output Feedback and Optimal 
Disturbance Suppression 


John R. Broussard 


CONTRACT NAS 1-1 7493 
MARCH 1986 


V1M10H1 A ‘U0I3WVW 

vsvn V.Hvaan 
^3IN30 H0HV3S38 A373NT1 



T? 


or Reference 


!,OT ;o EE tuoh xhis noor* 



IVI/NSA 




NASA Contractor Report 3968 


Extensions to PIFCGT: Multirate 
Output Feedback and Optimal 
Disturbance Suppression 


John R. Broussard 

Information & Control Systems , Incorporated 
Hampton, Virginia 


Prepared by 

Langley Research Center 
under Contract NAS1-17493 


NASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
information Branch 


1986 




FOREWORD 


The investigation described in this report was performed by Information 
& Control Systems, Incorporated (ICS) under Contract Number NASl-17493 for 
the National Aeronautics and Space Administration, Langley Research Center, 
Hampton, Virginia. This work was sponsored by the Applied Controls Branch 
of the Flight Control Systems Division. Mr. Richard M. Hueschen served as 
Technical Representative monitoring this contract. 
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ABSTRACT 


This report contains the theoretical development and the practical 
application of new control synthesis procedures for digital flight control 
systems. The new theoretical developments are the solution to the optimal 
multi-rate sensor output feedback problem and the solution to the problem 
of optimal disturbance suppression in the presence of windshear. Control 
synthesis is accomplished using a linear quadratic cost function, the com- 
mand generator tracker for trajectory following and the proportional-integral- 
filter control structure for practical implementation. Extensions are made 
to the optimal output feedback algorithm for computing feedback gains so 
that the multi-rate and optimal disturbance control designs are computed and 
compared for the Advanced Transport Operating System (ATOPS) . The perfor- 
mance of the designs is demonstrated using closed-loop poles, frequency 
domain multi-input sigma and eigenvalue plots and detailed nonlinear 6-DOF 
aircraft simulations in the terminal area in the presence of windshear. 
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1 . INTRODUCTION 


The proportional-integral-filter (PIF) control system is an established 
control structure which can be used to design and implement digital flight control 
systems for aircraft. The theoretical developments for the PIF design using an 
infinite time quadratic cost function, one sample rate and full state feedback is 
presented in Refs. 1-3. Successful flight tests of a PIF-like structure in 
position form is presented in Ref. 4. 

The PIF control system in a state space representation integrates well with 
the command generator tracker (CGT). The CGT, discussed in Refs. 5, & and 7 is 
a method for computing the trajectory plant states and controls follow when the 
output of a plant is tracking the output of a model. The command model can be 
used to generate coordinated turn paths for heading changes, smooth vertical 
height transfer paths for height changes or a complete three dimensional tra- 
jectory for curved path tracking. 

At present, all flight tested PIF control laws have been designed using 
full state feedback and the Linear Quadratic Regulator (LQR) approach. Problems 
have occurred due to the inability to accurately represent the actual aircraft 
dynamics and flight control system in the design aircraft model. For example, 
aircraft actuator dynamics, disturbance dynamics (gusts, wind shear), comple- 
mentary filter states, and analog prefilter states have not been included in 
PIF design models used to determine the full state feedback control law gains. 

The LQR approach would require feedback of all these states as well as the 
orginal aircraft states in an attempt to favorably, but impractically, alter 
filter and actuator dynamics as well as the aircraft's response to disturbances. 
The LQR approach guarentees + 60 DEG of phase margin and -6dB to + 00 dB gain 
margin, but may require feedback loops that are not practical to implement. 



The difficulties with full state feedback have been partially alleviated by 
extensive simulation of the PIF control system using a realistic truth model and 
iterative adjustment of quadratic weights in a more simple synthesis model. 

A fundamentally better approach is to use limited state or output feed- 
back in the PIF synthesis procedure. The realistic truth model can be used as 
the design synthesis model. Only practical measurements and feedback paths are 
used to design the PIF control gains. 

Until recently, the main obstacle to the use of the limited state feed- 
back approach has been the unavailability of a fast, reliable algorithm to 
compute the output feedback gains. This obstacle was recently removed by 
an algorithm whose derivation is presented in Ref. 8. This algorithm provides 
a fast, efficient and reliable method to obtain optimal output feedback gains 
for large order systems. 

Some progress in the use of output feedback, PIF and the CGT has been made 
in Ref. 3. Chapter 2 in the report begins with a more complete derivation of 
PIF with output feedback. The discussion in Chapter 2 provides the basic 
developments that the rest of this report will further investigate and genera- 
lize. 

The first generalization, discussed in Section 2A, is the use of the 
hierarchical or nested approach to design the feedback gains using optimal 
optimal output feedback. In the nested approach, one control loop can be 
designed with specific feedback paths, the designed control loop is closed, 
then the next control loop can be designed with a different configuration of 
feedback paths. Only a few small changes are required to convert the output 
feedback PIF design procedure to use the nested approach. The nested approach 
is a form of decentralized control, Ref. 9. 
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The theory of digital control design using PIF/CGT currently includes only 
the case where control effectors and sensor are all operating or sampled at the 
same rate. Due to practical considerations, sensor measurements, such as INS 
(Inertial Navigation System) signals, and external position measurements such 
as MLS (Microwave Landing System), GPS (Global Positioning System) and ILS 
(Instrument Landing System) may be available only at sample rates slower than 
the rates of onboard sensors used to stabilize inner loops. The space shuttle 
is an example of a vehicle which had to contend with sensor measurements avail- 
able at different rates. An important problem is to extend the PIF/CGT output 
feedback approach to the multi-rate sensor measurement case. The multi-rate 
sensor approach is particularly useful in designing inner-loop/outer-loop sensor 
control systems. The optimal multi-rate limited state feedback approach is 
derived in Chapter 3 and represents an important new contribution to multi-rate 
control synthesis. 

A contribution to the theory of control system design using the CGT was 
made in Ref. 5 in the research area of model following, and in the research 
area of disturbance suppression. Thus far, no flight tested PIF/CGT con- 
trol systems have exploited the disturbance suppression aspects of the CGT 
theory. 

An aircraft disturbance that has received considerable attention recently 
is windshear. In some types of windshear in the aircraft landing approach, 
the aircraft encounters a head wind, then a strong vertical downdraft, then 
a tail wind. Windshear is a type of disturbance that a flight control system 
must safely accommodate. 

A number of questions arise, however, that must be resolved in order to 
design a PIF/CGT control system with windshear disturbance assommodation. CGT 
theory assumes that the plant dynamics do not affect the disturbance dynamics. 
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Windshear disturbance models are different, in that certain aircraft states 
(height) affect the windshear dynamics. The first problem to be resolved then 
is to rederive CGT theory with the aircraft states coupled to the disturbance 
states. This derivation of a more advanced version of CGT theory is presented 
in Section 4B. A new type of matrix algebraic equation occurs during the develop- 
ment of the advanced CGT and is recognized as a generalized matrix Riccati equa- 
tion. A new, globally stable algorithm which finds a solution to the general- 
ized matrix Riccati equation is presented in Appendix A. 

Another CGT problem concerns the fact that the theory solves the distur- 
bance problem for perfect accommodation. Perfect or ideal disturbance accommoda- 
tion is accommplished by feeding forward the disturbance states in the control 
system. The disturbance states are not usually all exactly measured and the 
disturbance states not measured are sometimes difficult to estimate; A third 
question regarding the use of CGT theory concerns the use of imperfect, but 
adequate and practical, disturbance suppression. Given a PIF/CGT control 
system, which does not perfectly accommodate windshear, what is the effect of 
windshear on the aircraft states and how does changing the control gains alter 
the aircraft closed-loop windshear response? This question is resolved in Section 
4C. 

The steady state effect of stochastic disturbances on the closed-loop 
system is solved in Section 4C using a new concept called the stochastic star 
trajectory. The stochastic star trajectory shows that the steady state re- 
ponse of a plant driven by deterministic and stochastic disturbances decomposes 
into a deterministic plant response and a random or stochastic plant response. 

In Section 4D, a new type of optimal output feedback cost function is created 
using the stochastic star trajectory. The , quadratic cost is used to minimize 
the zero mean stochastic component and an algebraic matrix cost is added to 
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the quadratic cost in order to minimize the deterministic component of the 
steady state plant response. 

The last section in the report applies the new theories of multi-rate 
output feedback and disturbance accommodation to the design of an ATOPS air- 
craft path tracking autopilot. Nonlinear six degree of freedom simulations 
are made with the aircraft flying through wind shears. The energy probe sen- 
sor, Refs. 10 and 11, is included in the measurement vector and the effect 
on the design is presented. 
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2. OPTIMAL LIMITED STATE FEEDBACK PROPORTIONAL-INTEGRAL- 
FILTER/COMMAND GENERATOR TRACKER 


A derivation of the discrete-time output feedback proportional-integral- 
filter (PIF) control system combined with the command generator tracker (CGT) 
is presented in this chapter. A full state feedback derivation of the PIFCGT 
control system is presented in Ref. 2. The first derivation uses the usual 
integrated optimal output feedback approach. The nested approach, which takes 
advantage of the structure flexibility offered by output feedback is presented 
in the second section. In the nested approach, each control loop is designed 
individually in a sequential order. Each designed control loop is closed before 
the feedback gains for the next controller is computed. 

A. Integrated Approach 

The perturbation state vector, Ax, of the aircraft dynamics driven by control 
inputs, Au, and white Gaussian noise Aw is augmented to contain the perturbation 
control driven by the control rate, Av. Integral states A S are augmented to the 
state vector to operate on the aircraft output. Ay, 


Ax 


A 

0 

B_. 


Ax 


0 


Aw 

AS 

= 

H 

0 

D 


AS 

+ 

0 

Av + 

0 

Au 


0 

0 

0 


Au 


I 


G 


The control, Av, is used to optimize the quadratic cost function 


J = J o [Ax A£ Au] T Q 


Ax 

AS 

Au 


+ Av RAvdt 


( 2 ) 


The control rate, Av, is assumed to be constant over the sample interval. 

At. The continuous-time optimization problem is converted to the equivalent 
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discrete optimization problem using the sampled data regulator. Ref. 12. 

After converting the system dynamics and cost function, further simplifying 
assumptions and objectives are introduced into the optimization problem. Both 
the control position, Au, and control rate, Av, are required to be constant over 
the sample period. The discrete integrator A£^ is to be implemented digitally 
using Euler integration. The measurement noise present in the integration of 
Ay by the integrator is neglected. The objective of the control, Av, is changed 
from driving aircraft states to zero, to driving Ay, used in the integrator, to 
track the output, Ay^, of the command model: 

Ax . . . = $ Ax . + T Au , , (3) 

mk+1 m mk m mk+1 

Ay , = H Ax . + D Au , . (4) 

'ink m mk m mk+1 


In the derivation to follow, Au is assumed to change once at t and remain 

m o 

constant thereafter. In implementation, the command model control input is not 

constant and the command model dynamics are nonlinear. When Au^ is changing. 

Ay and Ay m become mismatched and their error is governed by the closed-loop 

dynamics and the integrator. When Au is constant. Ay eventually tracks Ay^ 

& & 

along the star trajectory , Ax and Au (assuming no plant parameter variations). 
The star trajectory is discussed in Refs. 1 to 3. The star trajectory is a 
linear system version of the nominal trajectory which is used when discuss- 
ing tracking for nonlinear systems . 

The star trajectory for discrete-time systems with a constant command model 
input is determined from 


f A 4 


A 11 

A 12 


' ^ 

E 

* 

< 

1 

1 

l> 

c 

s*r * 

ii 

A 21 

A 22 


Au 

L m J 


7 



( 6 ) 


The feedforward matrices A. . satisfy 


R*-d 

l 


pH 

pH 

< 

i 

A 12 


' A u ( V J) 

A ,, r 
11 m 

i — 

« 

D 


A 21 

A 22 


i 

D 

m 


which is a solvable matrix algebraic equation. The star trajectory is a con- 
venient notational abstract and is not generated in implementation. 

The tracking objective of the control law is introduced into the design 
by defining the variables 


Ax, 


~ 7 

= Ax^ - Ax£ ; Au = Au^ - Au^ 


Av, = Av, - A 01 (Ax . . . -Ax , )/At 
k k 21 m,k+l m,k 



and the discrete cost function. 


J ' k=-l 


; Ax£ QAx k + 2Ax k MAv k + Av £ RAv k 


(7a, b) 
( 8 ) 
(9a, b) 


( 10 ) 


At this stage, the star trajectory for the integrator is undefined, but is 
chosen as part of the optimization process. 

The star trajectory, by definition, must satisfy the plant dynamics 
with noise sources set to zero. Subtracting the star trajectory dynamics from 
the plant dynamics and using Eqs. 7 to 9 yields 


iS k+l" 


$ 

0 

r ‘ 


Ax, 

k 


O' 


" w k 


= 

AtH 

I 

AtD 


A ^k 

+ 

0 

Av, + 
k 

0 

- A Vi- 


. 0 

0 

I . 


.Au, 

k 


,‘IAt. 


.0 . 


$ r w, 

■ k 


( 11 ) 
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The plant is assumed to be tracking the model for constant Au^ previous to t =0. 

* 

The cost function starts at -1 since Au , and A? are to be determined. The 
A * 

quantities, Au and Ax are defined in Eq. 10 using Au instead of Au 

_1 _1 m,o m,-l 

The standard use of the linear quadratic cost function is to regulate nonzero 

initial condition states to zero. The cost function in Eq. 10 is constructed 

so that the non-zero initial conditions in Ax ^ are caused by the change in the 

command input at T- o = 0* The control law optimally transfers the system between 

star trajectories for a step change in u . In practice, the control law performs 

m 

well if u^ changes intermittently or "slowly" varies. 

To clarify the construction of the cost function, a brief scenario of a 

step change in Au^ is presented. At k=-2, the plant is tracking the star 

trajectory and the error quantities Ax and Au ^ are zero. At k=-l, the star 
* * 

trajectory for Ax and Au ^ is switched to the trajectory being generated by 

the new Au which will occur at k=0« The command model is chosen so that the 
m 

command model states at time k=0 are immediately affected by the command model 
control at k=0 as shown in Eqs. 3 and 4. The new star trajectory causes large 
errors to appear in Ax which is being weighted in the cost function. The 
control increment, Av , a variable to be optimally chosen, feeds back the large 
errors to generate Au q . The control Au q , generated from Au and Av ^ , immediately 
moves in a direction to reduce the tracking error when Au^ is changed at k=0. 

What this means in implementation, as shown in Eq^ -3 and Eq. 40, is that u m ^ is 
fed forward and directly affects u^. Similarly, starting the cost function at 
k=-l and computing the star trajectory at k=-l, using Au^ at k=0, aids in 
choosing A£_^ as discussed in the next paragraph. If the cost function were 
to start at k=0, (the standard optimal control starting point for the non-zero 


initial condition problems), then Av^ would be optimally chosen to affect Au^ . 
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The proper optimal value for Au^ and A£ q would be unclear; Au^ and A£ would 
become unknown initial conditions in the plant state vector. 

The next assumption used is that not all aircraft states are available 
for feedback. The states that are available for feedback are represented as 


Ay MEAS 


"c 0 

o” 


1 

> 

< 

i 


1 * 

> 

i* 3 

Au 

= 

0 I 

0 


Au 

+ 

0 

> 

vHti 

1 

k 

.0 0 

I. 


-A5- 

k 

. 0 _ 


( 12 ) 


The states observed from the aircraft model are corrupted by white Gaussian 
measurement noise Av^ with covariance, V. The control internal states Au 
and _A§ are noise free. 

The globally optimum solution for the problem construced thus far is a 
singular Kalman filter combined with the full state feedback gain from the 
linear quadratic regulator, Ref. 13. The globally optimum solution is not 
necessarily robust. Ref. 14, or straightforward to implement. Greater flexi- 
bility in the control design process is obtained if the quadratic cost is 
minimized using a prespecified control structure. The class of control laws 
considered in this chapter are restricted to be of the form 


A \ = - < K y K u V 


Ay, 

aG 

a£ 


MEAS 


= - K y. 


(13) 
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For the plant dynamics shown in Eq. 11, the feedback gain constraint 
shown in Eq. 13 and the following conditions: 


E l A ”ki ■ 0 E | i “ k4 "ki * '% 

E Ki - 0 e K % T i - vs kj 

EjAx Ax T I = X 
l _1 - 1 ) o 

e K a ^I ■ e K - e S a \ ■ 0 


the cost function in Eq. 10 is modified as follows (Ref. 15) *: 


J(K) - T W t + V 


J t ‘ J-l AS tk * Ax tk + 2 A ^k fi A "tk + A5 tk E A "tk 


J s ■ tj+T E J-l A ^k 5 Ax sk + 2 ^sk ** A; sk + A? sk E Aj sk 


J(K)= j tr J P(W + X Q )j + j trjK T (T T PT + fi.) K V 


The notation "tr" denotes the trace of a matrix. is the transient cost with 

noise sources set to zero while is the average stochastic cost. The trade- 
off between J t and J g is accomplished by varying X q with respect to W. The 
effect of X q is like adding pseudo process noise to the plant to improve con- 
trol system robustness. The benefits of using Eq. 18 as the cost function are 
discussed further in Ref. 15. 

The matrix, P, in Eq. 21, satisfies the Riccati-like equation 


— T — — T — T ** — — _ _ _ r r _ r r ^ t* 

P = $ CL P $ CL + C K R K C + Q - M K C - C K M 
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<I> is the stable closed-loop plant matrix 
L»L* 


= * - r K C (23) 

CiJL 

The necessary conditions for J(K) to have a minimum, are derived in Ref. 8 
The necessary conditions are: 

o There must exist a gain K so that is stable. 


o The gain K must satisfy 

(? T p T + r) k ( c s c T + v) = (r T p $ + m t ) sc t 


where 


(24) 


s = $> T s®?L + (w + x ) + r k v k t r T 

LL L<L O 


(25) 


A gain which satisfies the necessary conditions is not necessarily unique. 

A A 

The next objective is to determine A£ . The choice of A£ has a strong 

effect on Au , the first control effort that occurs when it is learned the 
o 

Au^ has changed. From a purely tracking viewpoint, the noise sources cause 
the states and controls to vary with zero mean about a trajectory which trans- 
fers the states and controls between star trajectories. The cost minimized to 
determine the state and control trajectory which intercepts the star trajectory 
is given by Eq. 19 where the noise sources are zero. Substituting the control 
system into the cost, J t , yields. 


= 




Ax^ k [Q - M K C 


Q 

_t _t _ r r «^ r r 

C K M + C K R K C] Ax 


tk 


(26) 
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Only the initial value of the state vector is unknown in the transient cost, 

CO 

k+l _ , \ k+1 


j -J,A? /?) sfi r , 

t k=-l t,-l V CL / x \ CL/ t,-l 


( 27 ) 


The infinite sum converges to a matrix, P, if $ is stable. The matrix P 

LL 

is the same as the matrix P shown in Eq. 22, hence; 


J t = Ax^ P Ax_^ = Ax^ 


P P P r 

XX xu xc, 


T 

P P P 
xu uu u^ 


„T „T 

M P u? 


Ax 


-1 


(28) 


For tracking without noise, Ax ^16 given by 


Ax 


-1 


Ax “ 


r * “1 

Ax , 

-1 


-1 

Au -i 

- 

A * 

AU -1 




5i_ 

” 

r-J 


(29) 


but 


p-i r 


' A ll 

A 12 _ 


7 

i 

i— 

i 

1 


_ A 21 

A 22. 


- A vJ 


(30) 


hence 


Ax 


-1 


A ( u , - u ) 
12 m-1 mo 

A 00 ( u . - u ) 
22 m-1 mo 


L < - < 


(31) 
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Treating the increment in the integrator as an unknown, £ = A£_£ - A£ , the 
value for £ that minimizes J'is given by 

3j/3£ = 0 (32) 

The solution is 


P 22 1 + (P H A 12 + P u{ V <A Vl - iu mo ) = 0 


(33) 


From Eq. 33, the value for A£ must be given by 


A£,, 


-[° a] 


Ax 


’mk 


(34) 


Au 


mk+1 


where 


A ■ - P 5C [ P *I A 12 + P u l A 2 2 ] (35 > 

Evaluating Eq. 34 at k=-l and k=0 and substituting into Eq. 33 demonstrates 
that the solution is correct. Using the feedback gain in Eq. 13, the per- 
turbation control system is 

Au^ = Au^_^ + At Av^_^ (36) 

Av k-1 (I " At K tP Av k-2 " K y (Ay k-l " Ay k-2 ) “ K £ (A ^k-l " A ^k-2* (37) 

The previous value, Av^_ 2 » is subtracted from Av^_j to obtain the incremental 
expression in Eq. 37. The perturbation variables and trim variables are elim- 
inated from the incremental expression using a large number of substitutions 
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and cancellations (discussed in Ref. 2). The use of output feedback in Eq. 37 
does not significantly alter the derivation in Ref. 2, which is for full state 
feedback. The implementable equations for the PIFCGT control law become, 


u, = u, . + At v. . + A„ , (x . - x , . ) 

k k-1 k-1 21 m,k m,k-l 


e k-l y MEAS,k-l C A 11 X m,k-1 


V k-1 - (! - At KJ v k _ 2 - K y (e k _ x - e k _ 2 ) - At - y^.,) 


(38) 

(39) 


K 


um 


+ (K r A + K C A. . + K A 0 „)(u m - u .) 

E, y 12 u 22 m,k m,k-l 


(40) 


The gain, K^, which feeds forward the command generator forcing function 
increment, is a linear combination of feedback and feedforward gains. The 
effect of K is to improve the transient response of the outer-loop control 
system by changing closed-loop system zeroes without affecting closed-loop 
poles. 

B. Nested Design Approach 

The output feedback control synthesis approach discussed in the previous 
section is an integrated control design. Each measurement is fed back to every 
controller. The feedback gains from measurements to controllers are all computed 
simultaneously. If a control actuator fails or a control actuator saturates, 
and becomes inoperative for a length of time, the control feedback network 
undergoes a structural change. The feedback gains designed simultaneously do 
not guarantee that the reduced controller situation will remain stable. 


Situations also occur where it is not desirable to feedback every measure- 
ment to every controller. In a jet transport, for example, feeding back pitch 
rate, normal acceleration and the height integrator output to the elevator, but 
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not to the throttle, can be a reasonable design constraint. 

The nested design approach allows the designer to specify the measure- 
ment vector for each controller. In addition, each control loop is designed 
in a sequential hierarchical manner with critical control loops designed first. 
The design of the nested approach begins by choosing a controller and the cor- 
responding desired measurement vector. The optimal output feedback optimiza- 
tion problem is solved to obtain the feedback gain. The control loop is closed 
to create a new plant model. The next controller in the hierarchy is used for 
design with a new measurement vector and the process is repeated. The con- 
trollers can fail or saturate in a particular order from the outer most con- 
troller loop to the inner loop and the closed-loop system remains stable. 

Since the development of the nested approach, two alternative techniques 
have been developed that are considered to be more useful for designing decen- 
tralized measurement feedback loops and accommodating control surface failure. 
An integrated approach to the decentralized measurement feedback loop pro- 
blem is presented in Ref. 16. An integrated appraoch to the combined decen- 
tralized control and control surface failure problem is solved in Ref. 17. 

The nested approach applied to the PIFCGT design begins as in Section 
2A up to Eq. 21. Computational considerations require that the design should 
proceed from control loop to control loop with no change to matrix dimensions. 

— — . A A A A - 

The design matrices $, T, C, Q, M, R, W, V are used to form a second group 
of design matrices where one controller is used for design and there are 
other controllers that have not been designed yet. Controllers not designed 
yet and integrators not used for feedback yet, cause $ to have uncontrollable 
states with eigenvalues at 1. 

The following operations are performed on the second group of design 

— “ — AAA — 

matrices C^, Q^, ’ v p to create a stabilizable optimization 

problem: 
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r - null columns in T for controls not being used 

C ^ - null rows in C for measurements not being used. These measurements 
include controls and integrators not desired. 

A A 

- not changed from Q 

A A 

- not changed from M 

A 

- not changed from ‘R "but must have full rank 

^1 - Diagonal elements of $ for control states not previously designed; 
integrator states not being used, are artifically stabilized by 
placing a small negative number in the diagonal position. 

- Positive scalars are added to diagonal elements of W for the same 
control and integrator states stabilized in 

— Positive scalars are added to diagonal elements of V for measure- 
ments whose rows are nulled in 

The changes cause to be stabilizable; control states and integrator 
states which are neutrally stable and not used in the design are made unobservable, 
uncontrollable and stable. The states that are uncontrollable and unobservable 

— A 

do not cause P. and S to reduce in rank. The changes to and R (if required) 

— /v — _x 

cause T^Pr^ + R and CSC + to be invertible so that Eq. 24 has a solution. 

None of the changes adversly affect the design of the control gains for the 
nested loop being synthesized. An alternative, but computationally unattrac- 
tive approach, would be to reduce the dimensions of the design matrices to 
eliminate unobservable and uncontrollable states. 

After a design for one control loop is completed the system matrices are 
altered as follows : 

5 = 5 - f x k x c 1 (41) 

A A rp /\ A T 1 T ^T 1 

Q = Q + Cj Kj R K 1 -C 1 - M Kj Cj - Kj M (42) 
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( 43 ) 


A A rp [| fTi 

w^ + ^^vKjrJ 

The alterations reflect the fact the is now a known matrix and must be 
absorbed into the design. Using the new system matrices, $ T c , ... etc. 
can be computed for the second controller in the hierarchy. The operations 
continue until all control loops are designed. When the nested design is 
completed by solving m different optimal output feedback problems, the total 
control system is 

v, = Tk, C, + K„ C n + ... K C x, = K C x, 
k 1 1 2 2 mmjk k 

The PIFCGT design requires that the P matrix defined 

* — — 
find A E, . The P matrix in Eq. 28 is equal to the P^ 

the last nested design performed to obtain K . 


(44) 

in Eq. 28 be used to 
cost matrix obtained in 
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3. OPTIMAL MULTI-RATE LIMITED STATE FEEDBACK 


A multi-rate control system has sensor measurements, computer calculations 
and control commands performed at a variety of different sampling intervals. 

The advantage of using a multi-rate control system is to better utilize the 
computational capability of the on-board computer. Critical control loops with 
fast dynamics can operate at fast sampling rates, while less critical control 
loops with slow dynamics can be processed at slow rates with little degradation 
in overall performance. Multi- rate control designs have been incorporated in a 
number of digital flight control systems for aircraft. 

Computational capacity is continually improving, making the need for multi- 
rate controls less attractive. A problem that will persist, however, is that 
all sensor measurements will propoably not be available at the fastest com- 
putational rate of the control computer. 

An aircraft, for example, using a full complement of on-board and external 
sensors, (such as the microwave landing systems (MLS)), often does not have 
these sensors all sampled at the same rate. Body mounted sensors are available, 
or may have the equipment potential to be available, at a fast sample rate, 
while INS and aircraft geographical position measurements are available at 
slower rates. 

A procedure that allows each sensor to be feed back at the sensor's 
sample rate is developed in this chapter using an optimal, multi-rate output 

J 

feedback synthesis approach. The increased computational capacity available 
when control commands are computed at different rates is also a feature of 
multi-rate sensor feedback. The control computations for slow rate sensors 
need only be performed at the sample rate of the sensor. The control calcula- 
tions for control commands to the actuators are assumed to be performed at the 
fastest control rate. Each sample rate of the different sensors is assumed to 
be an integer multiple of the fastest control rate. The assumption placed on 
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the sensors sample rate guarantees that the design model and cost function for 

the control design problem are periodic in time. 

Optimal multi-rate limited state feedback does not appear to have been 

previously investigated in the literature. The optimal full state feedback 

multi-rate control problem has been investigated in Ref. 18 using a state- 

space approach. The optimal control solution in Ref. 18 is a periodic sequence 

of full state feedback gains. The derivation in this report will take advantage 

of the structure flexibility offered by optimal limited state feedback and 

require that the multi-rate gain which minimizes the quadratic cost function 

be constant. A constant feedback gain reduces implementation complexity 

over a periodic sequence of feedback gains, particularly as the number of 

cycles in the period increases. 

A. Aircraft System Model and Cost Function 

The measurement vector is assumed to be separable into two groups, y^ and 

y . The y. measurements are sampled at the fast rate r_. The y vector is the 
s f f J s 

collection of sensors sampled at rates slower than the fast rate. The number 

of slower rates is assumed to be n. The n slow rates are grouped into the 

ordered sequence r , r „, ... r . The fast rate divided by the slowest 
si s z sn 

rate, r^/r , is Integer r. The sequence of fast and slow measurements 
repeats itself every r samples in a periodic fashion. A slow rate sensor output 
is held constant at the last sampled value during the fast sample times when the 
slow rate sensor is not being measured. 

The linear time- invariant model of the dynamics of the plant is represented 
at the fastest rate as follows: 

Vi ° * *k + r u k + “k (45) 
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The control is and x^ is the state vector of the plant model. The white, 
zero mean Gaussian process noise disturbing the plant is w^ with covariance W. 
The control system is assumed to have the fixed structure, 



(46) 


and K are, respectively, the fast rate constant feedback gain and the slow 
i s 

rate constant feedback gain. The measurement vectors are represented as 


II 

4-1 

C f \ + v f ,k 

k = 1, 2, 3 ... 

(47) 

i = 

1, 2, ..., n 



^si,k 

• C . x^ + v 
si tc si,k 

k = 1, r r /r . + 1 ,2 r./r . + 1 , ... 
f si f si 

(48) 

y si,k 

^si,k-l 

k ^ 1, r_/r . + 1 , 2 r c /r . + 1 , ... 
f si f si 

(49) 


The vector y g ^ is composed of the slow sensor measurements that are 

available at the slow rate r g ^* At the fast rate, Eq. 48 represents time 

samples when the y . slow rate sensors are measured. The y . measurement is 
si 'si 

held constant for all other fast sample times as shown in Eq. 49. The white 
Gaussian measurement noise sources are v^ ^ and v^ ^ with covarances and 


Equations 47 to 49 can be combined with the system model dynamics in 
Eq. 45 using the periodically time-varying matrix A^. The periodically time- 
varying matrix A^ is given by 
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0 


x 1 j 2 j * • « j n 


6 . , = I; k = 1, r./r . + 1 , 2 r /r . + 1, .., 

si,k f sx f sx 

0; k ^ 1, r r /r . + 1 , 2 r^/r . + 1, .. 

f sx f sx 


Using the periodic matrix, A^, the slow measurements can be combined into 
the measurement vector and represented at the fast rate as follows: 


y s,k+l \+l (C s *k+l + V s,k+1 


) + (i - A k+1 ) y g ^ k 


The combined plant and slow rate measurement periodic model has the follow- 
ing state space representation: 


O x 


\+l C S 1 y s. 


0 (I_A k+l ) J y s,k 


I 0 w, 


° Vl I * * * V s, 


The diagonal elements of A k switch between zero and one depending on 
whether the slow measurement should be represented by Eq. 48 or Eq. 49 at 

the kth sampling interval. The A k matrix repeats itself every r samples. One 

cycle of A k matrices are written as A^, A^, ... A^. The slow measurement 
observation matrix C g is composed of the matrices C g ^ to C gn packed rowwise, 
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( 54 ) 


C = 
s 




C 

sn 


Multiplying Eq. 53 by the matrix 


A. c 

k+1 s 


, the inverse of the 


left partition matrix in Eq. 53, changes the plant representation to a standard 
time-varying difference equation format:: 


k+1 


x 


k+1 

^s ,k+l 


A k+1 C s* ^k+l 


x k 


r , 

rj 

+ 

T 

n<" 

cn 


_ A k+! 


c r 


u k + 



(55) 


The feedback control system is represented as follows using the periodic 
model states: 



(56) 
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The cyclic nature of the matrices shown in Eq. 55, means that there are 
recurring moments in time where the matrix cycle begins and ends. Substituting 
Eq. 56 into Eq. 55, the periodic closed-loop plant representation for a cycle 
which begins at index k and ends at index k+r-1 is 


X k+1 

^+2 





KC) X k + E 1 \ - F 1 K \ 

KC) \+l + E 2 W k+1 ” r 2 Kv k+1 


KC) 


k+r-1 


+ ■ E w, 


r k+r-1 


- T 'Kv. 


k+r-1 


where 



$ 

0 


$ 

0 

*1 = 

A 0 C $ 
2 s 

!-A 2 

V# • 

■e 

N> 

II 

A.C $> 
3 s 

!- A 3 

" $ r-l = 

$ 

0 

; $ = 

$ 

0 

A C $ 
r s ' 

I-A 

r 

r 

A.C $ 
1 s 

1 

»— H 
<1 

i 

w 


(57) 

(58) 

(59) 


(60) 


The other matrices in Eqs. 57 to 59 are similiarly defined. 

The multi-rate optimal output feedback cost function can be constructed 
in a variety of ways. The following simple approach is used in this develop- 
ment. The sampled-data regulator cost function at the fast rate in standard 
form is given by 


= lim 

N-*» 



( 61 ) 
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The weighting matrices Q, M, and R can be obtained from a continuous quadratic 
cost and plant model using the fast sample rate and a suitable software package 
such as ORACLS, Ref. 19. The additional state, y g caused by the slow measure- 
ments, is arbitrarily introduced in the cost with no weight 



(62) 


(63) 


The y g states are a part of x and can be weighted using Q. The design objective 
is to determine the control system gain shown in Eq. 56 which minimizes the 
cost shown in Eq. 62. 

It can be shown that if a periodic system is stable, then the covariance 
of the states of the periodic system reach a periodic steady-state. Assuming 
the periodic system shown in Eq. 55 can be stabilized using Eq. 56, the peri- 
odic steady-state covariances are the solution to the following sequence of 
equations : 

5 2 = (^ - Tj K C) S 1 ($ x - r x K C) T + E x W eJ + r x K V K T r[ (64) 

5 3 = (< ^2 “ r 2 K C) s 2 ($ 2 - r 2 K C ) T + E 2 W E 2 + r 2 K v ^ V 2 (65) 

s = ($ - r k c) s ($ - r k c) T + e w e t + r k v k t r T (66) 

lrr rrr rrr r 
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The periodicity is shown in Eq. 66 where the covariance matrix equation 
cycle begins to repeat itself. In periodic steady-state, the cost function 
can be rewritten as 


N 


J = lim 


Nhw) 2rN 
s s 


£ ^cycle 2r ^cycle 


i=0 


( 67 ) 


The integer, N , is the number of periodic cycles in N. 
s 

the weighted states and controls for one cycle. 


J , is the sum of 
cycle 


cycle 


=- tr 


j (Q - M K C 


T T T T T T 

C K M + C K R K C) S x + S 2 + . . . + S 


J! 


4- r tr | K T R K V | 


( 68 ) 


B. Necessary Conditions For An Optimal Controller 

Adjoining the equality constraints shown in Eqs. 64 to 66 to Eq. 68 yields, 


^ 2r ^cycle + 

t» (-s 2 + (® - Tj k c) s x ($ : - r x k c) T + Ej w eJ + r L K V k t rj) P P + ... 

+ tr (-s. + ($ - r k c) s ($ - r k c) T + e w e t + r k v k t r T ) p T (69) 

1 rr rrr rrr rr 

The Lagrange multipliers are P^ to P^. The conditions necessary for J to have 
a minimum are aj/ 8S x = 0 ... 3J/3S r = 0, = 0 ... 3J/3P r = 0, 3J/3K = 0. 

The 3J/3S^ sequence determines the equations for P^: 

P r = ( ^1 " r i K C)T P 1 ($ 1 - r ! K C) + Q (70) 
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(71) 


P 1 = ( $ 2 ” r 2 K C )T P 2 ($ 2 ~ r 2 K C > + Q 

P r-1 = ( *r " F r K C)T P r ($ r " r r K c ) + Q (72) 

The necessary condition 9J/3K determines the equation for the feedback gain 
K; 

rk(cs, c T +...+cs c T ) + (r! p, r, + ... r T p r + r r) k v + 

1 r ill rrr 

r p p, r. k c s c T + . . . + r T p r c s c T = rT p. $ s c T + ... + 

111 1 rrrr 111 

T T T T T T I 

r P $ s C + M S. C + ... S c (73) 

rrrr 1 r 
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4 . NUMERICAL APPROACHES 

The numerical solution for and P^ given a feedback gain matrix, K, 
are obtained by expressing one of the covariance matrices and one of the 
cost matrices at the slow rate. Starting with Eq. 64, S£ can be substituted 
into Eq. 65 and so on until Eq. 66 is reached. When Eq. 66 is reached the 
only unknown is the matrix S^, 

S 1 = *1 S 1 ^1 + W * (74) 

where 


$ = ($ - r kc) ($ . - r . k c) ... ($. - r. k c) 

J r r r-1 r-1 j j 


(75) 


*r+l = 1 


(76) 


W. = E. W eT + r. K V K T rT (77) 

J J J 3 J 


* 

w 


r 


l 

3 = 1 



W. 

3 3+1 


(78) 


The covariance equation shown in Eq. 74 is easily solved. Once S^ is 
known, the other covariance matrices can be reconstructed from Eqs. 64 to 66. 

A 

The matrix is the representation of the plant at the slowest rate. The 

A 

periodic system is stable for a given K matrix if has eigenvalues within 
the unit circle. 

The solution for the cost matrices, P_^, are similarly determined : 


P = i T P $ + Q* (79) 

r r r r x 
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( 80 ) 


i T = ($ - r k c) T ($. - r k c) T ... ($. - r. k c) t 

3 1 1 2 2 3 3 


l l =1 


* V ~T A ~ 

q - A . q . 

J=1 J-l J-l 


(81) 


(82) 


The remaining numerical problem is the solution for K in Eq. 73. Equation 
73 is a Lyapunov-like equation with more than 2 entries, i.e.. 


A 1 X B 1 +A 2 XB 2 + A r+2 X B r+2 = C 


(83) 


The method currently being used for finding X is to use Kronecker products. 



(84) 


and solve for x by inverting the Kronecker product matrix sum. The Kronecker 
product matrix sum yields a symmetric, positive definite matrix if K is the 
optimal feedback gain. The vector x consists of the columns of the matrix X, 
stacked. Likewise, the vector c is a stacked version of the matrix C. 

A. Numerical Algorithm for Solving Necessary Conditions 

The proposed algorithm is a straightforward extension of the convergent 
algorithm discussed in Ref. 8. 

/\ 

1. Choose a starting gain K q which stabilizes Choose a positive 

scalar a < 1, so that J(K^) is decreasing in the iterative procedure. 

2. Solve for S^ in Eq. 74, then iteratively compute S£ through starting 
from . 

3. Solve in Eq. 79, then iteratively compute to P^ starting from 

P . 
r 
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4. Solve for K in Eq. 73. 

5. Find the next gain in the sequence 

K i+1 = K i + a ( K "V (85) 

6. If J (K^ + ^) - J(K^) and || 9J/9K^ || are less than some convergence 
criterion stop, otherwise set i = i+1 and repeat the sequence starting 
from 2. 

If the algorithm has convergence difficulties or the Kronecker product 
matrix is not positive definite, reduce a and start over from the most recent 
stabilizing gain as discussed in Ref. 8.. At present, conditions which guarantee 
algorithmic convergence for a>0 are unknown. 

B. Application To The PIFCGT Design Approach 

The multi-rate optimization problem is performed at the fast sample rate. 
Incorporating the command generator tracker into the design problem along with 
integrators and control difference weighting requires careful consideration. 

The derivation which follows results in the PIFCGT control system that can be 
implemented to control nonlinear dynamics and use multi-rate feedback. Assump- 
tions have been made to determine a solution. The integrator is computed at the 
fast rate but can use fast or slow measurements in the integration process. The 
control position is always updated at the fast rate. The command generator 
tracker and star trajectory are assumed to be computed at the fast rate. The 
measurement system does not affect the tracking objectives. A multi-rate com- 
mand generator tracker is possible, but is beyond the scope of this effort. 

Starting from the plant dynamics at the fast rate, the dynamic equation 
for the error variables x^ and u^ is determined exactly as in Eqs. 7a, 7b and the 
v^ dynamics are adjoined to the model: 
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( 86 ) 


VI 


$ 

f 


X 

+ 

~ 0 

“J 

k+1 

0 

I 


u 

k 

Atl 


The slow rate measurements tracking command error is defined as follows : 


sk 


y sk-l + \ (C 


s X k 


" y sk-l } 


( 87 ) 


The variable y , is defined by the above equation. The plant state, y , and 

SK S 

ft ^ ft 

star trajectory, y , cannot be separated for the y , variable because y 

S SK S 

is undefined. The state and star trajectory can be separated for C . 



C . x, = C . x, 
si k sx k 


C . 
sx 


X 


* 

k 


The integrator for slow and fast measurements is 

e k+l ■ + At H s 5 sk + At (y fk - y fk> 

* 

Subtracting t, from both sides yields 

5k+l ■ \ + At H s y sk + it H f ^ + it D f D k 

where 


( 89 ) 


( 90 ) 


y fk = H f \ + D f u k 


( 91 ) 


Grouping everything together, the multi- rate PIFCGT design model is 


X 



0 

T 

0 


X 


0 

1 


AtH f 

I 

AtD f 

AtH 

s 


l 

+ 

0 

u 


0 

0 

I 

0 


u 


Atl 

Js_ 

k+1 

\+l C s * 

0 

\ + i ° s r 

I_A k+l_ 


_ y s 

k 

0 
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I 

0 

0 


Vl C s 


0 

AtH 

0 

0 



(92) 


The cost function is constructed similiar to Eq. 62. The measurement equation 
is 


V 


~ c f 

0 

0 

0 


X 


V fk 

l 

= 

0 

I 

0 

0 


l 

+ 

0 

u 


0 

0 

I 

0 


u 


0 

- ? s- 


_0 

0 

0 

I_ 


-v 

1 

_0 


(93) 


Optimizing the cost 


function produces the constant gain feedback control system 


v. = K. K r K K 
-He L f S U s_ 


(94) 


The next step in the derivation is to increment the measurements. The deriva- 
tion for incrementing the measurements, which involves simple algebraic mani- 
pulations and careful attention to time indexes, becomes lengthy if more than 
one slow rate is assumed. The derivation continues with only two rates, r^ 

and r where 
s 

r /r. = r (95) 

s t 
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The generalization to more than one slow rate is straightforward. The slow 
measurement increment for y . using the definition in Eq. 87 is 


y sik y si 


, , = (C x. - C x. ) - (C x, - C x, );k=l, r+1 , 2r+l (96) 
,k-l s k s ic-r s k s k-r ’ ’ * 


y sik y si,k-l 


k/1, r+1, 2 r+1 (97) 


The error increment 

r~ 


e k " e k-l 


Y fk y f,k-l 


y sk y s,k-l 


(98) 


can be expanded to 


e. - e. , 
k k-1 


y fk ~ y f,k-.l 


— 

C A (x - x 

f 11 mk m,k-l 

\ (y sk ■ y s,k-d 


A, C A, . (x , - x ) 

k s 11 mk m,k-r 


C,. A (u - u ) 

f 12 m,k+l mk 

A, C A, _ (u . . , - u , , , 

k s 12 m,k+l m,k-r+l 


(99) 


The derivation for finding £ parallels the derivation in Section 2A up: to 
Eq. 33, ‘which becomes 


^-i L x? -i 


+ P r u + P_ y 
£u -1 s 




( 100 ) 


The star trajectory can be separated from x ^ and u ^ but not from y^ 

1? [ P *S (X -1 ' ’‘-d + p 5u (u -i - u -d + p 5y y s ,-i] <101) 


U - 5 -l ' - P C5 
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A ?s,-1 - ^3,-2 + A -1 C s (x -1 - x -l> 


( 102 ) 


Since only two rates are considered, A_^ is either 0 or I. If the matrices 
A and B are defined to be 


A ~ P CC (P x£ A 12 + P £u A 22 } 


(103) 


B = -P * (P- C A ) 

v, Cy s 12 


(104) 


then 


£ = £ + A (u -u . . ) +. A . B (u ■ -u ) 

-1 -2 mo m,-l -1 mo m,-l 


(105) 


If A j =1 then,. Eq. 102 is substituted into Eq. 101 to yield the B matrix in 
Eq. 104. If A_^ = 0 then, the term y g ^ should be used in Eq. 101, but 
y __ is zero since the control system is assumed to be tracking the star 
trajectory for k<-l. 

Using a time index, the integrator star trajectory is defined using 
Eq. 105 as 


? k ^k-1 (A + A k B) (u m,k+l" U mk ) 


(106) 


If u^ is changing then is periodic if the B matrix is nonzero. If u^ 
* 

is constant then £ is constant. 

K 

Incrementing the integrator and using the integrator star trajectory 
yields 


L - I. . = L - L , - A (u , . . - u ) - A, B (u , . - u , ) (107) 

k k-1 k k-1 m,k+l m k k m,k+l m k 
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where 


^k-1 At H s y s , k- 1 + At (y f,k-l y mf ,k-l^ 


(108) 


Defining the variables 


y £sk H s y sk 


(109) 


y m sk H s y mk 


( 110 ) 


then, from Eq. 87, 


5 5sk ' (I - V Wl + \ (y 5sk " 5 'msk ) 


( 111 ) 


where A^ is a matrix that satisfies the equation 


H A, = A H 
s k k s 


( 112 ) 


If A, is a zero matrix then A, is a zero matrix. If A, is a S> x £ 
k k k s s 

identity matrix then A is an mxm identity matrix. 

rC 

Grouping all the results together the multi-rate output feedback PIFCGT 
implementation equations are 

(from Eq. 86) 


u, = u, + At v, . + A. . (x , - x , , ) 
k k-1 k-1 21 m,k m,k-l 


(113) 


(from Eq. 99) 


6 k-l “ \-2 


y f,k-l y f,k-2 C f A 11 (x m,k-l " X m,k-2 ) 

A k-1 (y s,k-l “ y s,k-r-l ) “ A k-l C s A ll ^m.k-l^mik-r-P 


(114) 
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(from Eq. Ill) 


^s;,k-2 (I A k-2 ) y £s,k-3 + A k-2 (y £s,k-2 y ms,k-2 ) 

(from Eqs. 94 and 99) 

^k-1 = (I + At V — k-2 + [ K f K s] [-k-1 " %-2^ + 

At K ? (y £s,k-2 + (y ?f,k-2 “ y mf,k-2 )) + 

(-K ? A - ^ B - K f C f - K u A 22 ) (u*- u^) + 

" K s A k-1 C s A 12 (u mk r “m.k-r^ 

Although the feedback gain is constant and the CGT feedforward gains are 
constant, the implemented feedforward gain is periodic because of the in- 
tegrator star trajectory. 


(115) 


(116) 
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5. DISTURBANCE ACCOMMODATION FOR WINDSHEAR DYNAMICS 


One of the primary purposes of a flight control system is to stabilize 
the aircraft about a desired trajectory in an enviroment where external dis- 
turbances are affecting the aircraft dynamics and sensor measurements. Some 
of the disturbances are difficult to model and are best characterized by zero- 
mean Gaussian white noise. Other types of disturbances are best described by 
differential equations driven by Gaussian noise. Examples of the latter type 
of disturbance include the well known Dryden model for gusts, steady-state 
atmospheric wind and windshear. The objective of the control system is to 
stabilize the plant while having the capability to cope with the disturbances. 

Disturbance accommodation has a history of investigation with progress 
made by a number of researchers. The procedure for computing the control 
system in this paper is to minimize a linear quadratic cost function. Efforts 
in this area include Refs. 20 to 23. The objective in the references in to 
construct two signals in the control variable, u, i. e., 

Au = Au + Au (117) 

f a 

where, Au f , is designed to stabilize the plant while the other signal, Au , 
is chosen to accommodate the disturbance. One method for determining Au^ 
exactly using feedforward control is discussed in Refs. 5 and 24. The 
effect of the disturbance is eliminated on selected system outputs. Feed- 
ing forward the disturbance in Au usually requires some form of dynamic com- 
pensator (observer, servo compensator), since the disturbance is rarely mea- 
sured. The feedback system for Au^ usually does not contribute towards steady- 
state disturbance suppression. 
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An alternative to Eq. 117 is investigated in Refs. 25, 26, and 27 where 
only the Au^ signal is used to stabilize the plant and eliminate the effects 
of the disturbances on selected system outputs. The disturbances are suppressed 
using feedback gain eigenstructure assignment. 

An alternative to the approaches just discussed (which can be considered 
more practical) consists of constructing a feedback controller which strives 
to minimize, rather than completely eliminate, the effect of disturbances on 
selected system outputs. Investigations into the more practical approach have 
been made in Ref. 27 using pole assignment to achieve stability and Ref. 28. 

The disturbance accommodation procedures in Refs. 21 to 28 all assume that 
the disturbance dynamics are not affected by the plant dynamics. Models for air- 
craft windshear have a plant representation where the aircraft plant states affect 
the windshear dynamics as discussed in Ref. 29 and presented in the next section. 
Changes in aircraft height cause the windshear state to change in value. 

In the rest of this chapter a new tack is pursued for computing disturbance 
suppression control systems using the optimal control approach. The designer 
can specify any control structure that can be modeled using output or limited 
state feedback. The chosen control structure may or may not be able to exactly 
suppress the deterministic disturbance. The plant model and disturbance model 
are expressed in discrete time. The plant model and disturbance model can be 
completely coupled and each model can be driven by Gaussian white noise. The 
objective is to design the control system gains optimally so as to minimize 
the stochastic effects of the disturbances on the plant states and minimize 
the deterministic effect of the disturbance on selected plant outputs. The 
resulting cost function used to determine controller gains has two distinct 
parts. 
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The dual cost is accomplished using a result presented in Ref. 30 and devel- 
oped further in this chapter. The steady-state plant state and control response 
can be decomposed into a deterministic component dependent only on the distur- 
bance which satisfies a differential equation and a zero mean stochastic re- 
sponse due to the fcuassian white noise sources. The necessary conditions for 
the dual cost function to have a (local) minimum are derived in this chapter. 

A numerical algorithm is also developed to solve the optimal disturbance mini- 
mization problem necessary conditions. The special case where the disturbance 
is not affected by plant states is presented. A design example using PIFCGT 
is discussed in Chapter 5. 

A. The Windshear Model 

A basic wind model relationship in the aircraft's local level plane at 
a position above the earth's surface for windshear was developed by C. Belcastro 
and A. Ostroff at the NASA Langley Research Center, Ref. 29. The NASA Langley 
windshear model is presented in this section. The wind shear plus steady wind 
affecting the longitudinal dynamics of the aircraft is modeled as 


— — 
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Defining the following variables 



W 

w 


w 

— ce 


U 

wc 


W 

wc 

- e 


(119) 


( 120 ) 
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( 121 ) 


The relationship can be restated as 


w = w + s z 
e ce z e 


is the wind velocity vector in the earth-fixed reference frame, w^ is 
the constant wind velocity, s is the wind shear gradient vector and z is the 
vertical height of the position in the atmosphere above the earth's surface. 
The wind velocity in aircraft body axis is given by 


w, = H A w = H A w + H„ s z 
b 0 e 0 ce 0 z e 


where 


c9 -s0 
S0 C0 


and 0 is the aircraft Euler pitch angle. Taking the derivative of Eq. 123, 
yields. 


w = 9 L w + H w +0 L n s z + H a s z + H A s z 
b 0 ce q ce 0ze 0ze 0ze 


where 


-S0 

-C0 

C0 

-s9 


Equation 125 can be simplified by factoring 0 Lg and substituting Eq. 123, 


w = 9 I w, + H A s z + H Q w + ‘H a s z 
b 0b 8ze 0ce 0ze 
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where 


X 0 = 


0 -1 

1 0 


(128) 


The perturbation representation of Eq. 127 is 


Aw, =0 I. Aw, + I w A0 + L q s z A0 + z H q As 
b o 0 b 0 bo 0o zo eo eo 0o z 


+ H Q s Az + L q w A0 + H q Aw + z L q s A0 
0o zo e Uo ceo 0o ce e 0o -zo 


+ z Hq As + Hq s Az 
eo Uo z Uo zo e 


(129) 


The o subscript indicates the variable is representing the nominal value at 

the desired flight condition. The vectors w and s are assumed to be con- 

ce z 

stant, i.e. 


w = s = Aw = As =0.0 
ce z ce z 


(130) 


The windshear gradient perturbation vector is assumed to be 0, As =0 

z 

Substituting Eq. 130 into Eq. 129 and simplifying, results in the following, 


Aw, = 0 Iq Aw, + Iq w, A0 + z L q s A0 + H q s Az 
b 0 b 0 bo eo 0o zo 8o zo e 


(131) 


Eq. 131 can be rewritten as 


Aw = A Aw, + H Ax + H. Ax 
b w b x x 


(132) 


where 


A = 0 1 

w 0 


(133) 


(2x2) 
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(4x5) 


0 0 z L s 0 

eo 0o zo 


0 0 0 


0 0 0 I Q w, H q s 

0 bo 0o zo 


(4x5) 


0 0 0 


Ax^ = ^Au Aw Aq A0 Az g J 


Grouping the wind velocity model in Eq. 132 with the standard representation of 
the plant dynamics produces, 


I 0 Ax 


•H . I Aw 
x b 


A D Ax 

w 

H A Aw, 

x w 


The matrix D is obtained from the Au and Aw column vectors in the A matrix, 
w 


D = [-A -A 0 0 ; A = A A A A a A 

w L u w J Luwq0 zej 


Multiplying Eq. 137; by 


G, *J 


produces the final plant representation 


H +H.A A +H.D Aw, 
xx w x w b 


The wind vector, Aw^, is affected by the plant states but not by the control 
states since it is easily shown that ILB is a zero matrix. Representing the 
aircraft plus wind dynamics in discrete time, however, causes the discrete 


control effect matrix T to have entries in all locations. 
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B. Deterministic Disturbance Accommodation Using The Command Generator Tracker 
The objective of this section is to determine the plant star trajectory 
when the plant output, y, is tracking a model output, y^ and the plant dynamics 
are affected by deterministic disturbances. The combined plant model and dis- 
turbance model are represented as follows 


Ax 


$ 

D 

w 


"Ax| 

+ 

"r ' 

As 


H 

L W 

$ 

wj 


As 

1_ 

r 

L wJ 



(140) 


A >"k ■ H x Ax k + H s As k 


(141) 


The aircraft perturbation state vector is Ax , the perturbation control vector 

xC 

is Au^> the disturbance vector is As^ and the selected perturbation outputs 
in which the effect of the disturbance is to be suppressed is Ay^. The out- 
put of a command generator, Ay^, is generated by the command model 


Ax , , = $ Ax + F Au 
m,k+l m mk m m 


Ay , = H Ax , + D Au 
v mk m mk m m 


(142) 

(143) 


The command generator state is Ax^ and Au^ is the command generator con- 
trol that is assumed to be constant. The dimensions of Ay Ay and Au are 

• m 

assumed to be equal. The aircraft output tracks the model output while sup- 
pressing the disturbance along the star trajectory. Using the same structure 
as in Chapter 2, i. e. a linear relationship, the star trajectory is con- 
structed as 


> 1 
X 

h 





“i 

Ax mk 

+ 

V 


Au 

u m ^ 


- S 23- 


As, 


(144) 
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Equation 144 is the solution to the tracking problem if the S „ matrices 

exist. The unknown feedforward S.. matrices can be shown to satisfy solvable 

i 3 

* * 

algebraic relationships by expressing Ax^^ ” Ax, first using the plant 
dynamics, Eq. 140, then using Eq. 144 and equating the two expressions for 

k k 

Ax^ + i - Ax^ Withput going through the details, the algebraic relationships 

for the S.. matrices are as follows: 
ij 


(*-i-Si 3 h u ) (r-s 13 r v ) 


H 


S 11 S 12 


- S 21 S 22 


— 



Mb — 


$ 

r 


S 13 


H 

0 


_ S 23_ 



S,, $ - D 

13 w w 


-H 


+ 


S., H 
13 w 


11 m 


H 

m 


S.„ r 
13 w 


S, , 

r 


11 

m 


D 



m 



- 

1 


S 13 


S 23 


(145) 


(146) 


A numerical solution for and S 22 in Eq- 146 is obtained by inverting 
the quad partition matrix 



r" 

2 . 

' a n 

a i 2 " 

H 

o_ 


1 — < 

CM 

n 22- 


(147) 


and multiplying Eq. 146 by the inverse: 


13 


23 


"ll ° 


Q 21 ° 


'13 


23 


$ + 
w 


n il ° 


^21 ° 


'13 


23 


k r »] 


13 


23 


Q..D + O H 
11 w 12 s 


-°2l D v + n 22 H sJ 


(148) 
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Equation 148 is a generalized Riccati equation. A convergent algorithm for 
solving Eq. 148 has been developed and is described in Appendix A. There may 
be more than one solution or there may be no solution to Eq. 148. After 
and are determined, they can be substituted into Eq. 145 and the result 
can be numerically solved using the technique discussed in Ref. 6. The con- 
trol tracking problem solution in Eq. 145 is not independent of the distur- 
bance rejection solution because of the coupling matrices and r . 

Any control system of the form 

Au k = Au k +KC x (Ax k - Ax k ) (149) 

can be ideally used to stabilize $ and eliminate the steady-state disturbance 
in Ay, . In the control system, K is the feedback gain and C is the observa- 
tion matrix for the system states. Equation 149 can 'be rearranged as 
-f -a 

A„ k - K C x Ax k ♦ (S 23 -K C x S 13 ) As k ♦ <S 21 ~K CM^) Ax^ + (S^-K Au^ (150) 

which is the form shown in Eq. 117. If ^ 22 “^ ^ as of full rank, the distur- 

bance state. As, , must be measured or observed in order to implement Eq. 150. 

It is evident from Eq. 150, that if K is to stabilize the plant and eliminate 
the disturbance in the Ay, response without using a u control term, then K 

K ‘£L ' 

must satisfy the equation 

S 0 „ - K C S = 0 (151) 

x 13 

A feedback gain which satisfies Eq. 151 does not always exist. Necessary and 

sufficient conditions for the existence of a gain K which satisfies Eq. 151 are 

given in Ref. 26 for the case where H and V are zero matrices. 

0 w w 

The next section addresses the question: What is the steady-state plant 
response to the disturbance and white noise sources if the control system shown 
in Eq. 150 is implemented? The formulation includes the case of an estimator 
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observing As^, since the estimator dynamics can be adjoined to the aircraft 
dynamics to extend the state dimension of Ax. 

C. The Stochastic Star Trajectory; 

The command generator tracker and the star trajectory have always been 
applied to deterministic models. In this section, noise processes are intro- 
duced into the models and the concept of the star trajectory is used to in- 
vestigate the tracking response of the stochastic system in steady-state. The 
objective to determine a design technique for accommodating windshear dynamics 
that are driven by random noise. 

The aircraft model and disturbance dynamics shown in the previous section 
are generalized further and represented as follows: 

Discrete Perturbation Aircraft Model: 

Ax, - = $Ax, + TAu, + D As, + w, (152) 

— k+1 — k ^-k w — k — k 

Disturbance: 


As, . , = $ As, + H Ax, + r Au, + o, 
—k+1 w — k w — k w — k — k 


(153) 


Measurement : 


Az, = C Ax, + C As, + v, 
— k x — k w — k — k 


(154) 


Control System: 

Au, = - K Az, = -KC Ax, - K C As, - K Av, 

tC K. X — “K W — K. *"”iC 


(155) 


The zero-mean, white Gaussian noise sources in the above models have covari- 
ances:; 


E hc4 = w 


E 2k <Zk “ w s 


E 5k 2k =W sx 


(156) 

(157) 

(158) 


E 


= V 


(159) 
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All other cross covariances between the disturbance and measurement noise sources 


are assumed to be zero. The C matrix in Eq. 155 allows for the fact that some 

w n 

sensors, such as an airspeed sensor, measure a combination of the aircraft and 
disturbance states. The disturbance and plant process noises are allowed to be 
cross correlated as shown in Eq. 158. 

Given the control system structure shown in Eq. 155, the steady-state 


response for Ax^ is assumed to satisfy the equation 




(160) 


The assumption used in Eq. 160, where an infinite sequence of unknown feed- 
forward matrices is used to describe the star trajectory, has been used before 

* 

in Ref. 5 . In the final answer for Ax^» t * ie infinite sequences will be re- 
placed by a difference equation once the A^, B^, and matrices are known. 
The solutions for S^, A's, B's, and C's are determined by computing 

* A 

Ax^ + ^ - Ax^ as discussed in the previous section. Using the plant dynamics, 

* & 

the expression for Ax^^ - Ax^ becomes 


Ax* + i - Ax* = ($-r K C x - I) Ax* + (D w -T K C w ) As fc -r K Av fc + Aw fc (161) 

^ & 

Using 153, 155, and 160 Ax . — Ax , results in 

-k -k 
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- Ax* = S„. [($ - r K C 

k+1 k 31 L w w 


+ S 14 a k + 


! v> - X ] 

[" A l V A 2 A 2“ A 3 ‘"I 


As, + S,. (H 
k 31 w 


r K C ) Ax. - S. , r K Av. 
w x' k 14 w k 


Aw, 


Aw, 


k+1 


Aw, 


k+2 


- B 1 b t b 2 b 2 -e 3 . 


Av, 

+ [ G 1 Cj> c 2 c 2 c 3 ...] 

— » ““ 

A a. 

-k 

k 

A \+i 


A Vi 

A \+2 


A °k+2 

• 


___ « — 


(162) 


Substituting the expression tor Ax^ from Eq. 160 into Eqs. 161 and 162 and equating 
the two expressions yields solvable equations for the S^, A^, IL, and CX matrices. 
Defining the matrices 


4> - (t-TKC) -S„ (H - T K C ) 

c x 7 31 w w x' 


x = TK - S_. r K 
A c 31 w 


(163) 

(164) 


expressing Eq. 160 in the z-domain and substituting the solutions for the A_^, B ^, 


and matrices into Eq. 160 yields: 


Ax = S As , . 

31 l c c c 


[ -1-2 -3 2 1 

$ + $ z + $ z + . . . 

. c c c -* 

+ k l 


Aw 


X„ + X„ z + C 3 X„ z 2 + • • • _ 


c c c 


c c 


Av 


S„. + .$ 2 S„. z + <T 3 S,. z 2 + 

31 c 31 c 31 


-] 


A a 


(165) 


satisfies the generalized Riccati equation 
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($ - TK C ) S 01 + s 01 ($ - r K C ) - S 01 (H - r K C ) S 01 + (D -r K C ) =0 

x 31 31 w w w 31 w w x 31 w w 

Equation 166 can be expressed in the form of Eq. 200 of Appendix A if (0 - T K C ) 

X 


is invertible. 

The infinite z-domain sequences in Eq.165 have the closed-form relationship 


[i - z *J l - [i - 1; 1 * * z 2 + ...] 


(167) 


Substituting Eq. 167 into Eq. 165 and converting back to the time domain yields 


Ax k = S 31 As fc + Ax k (168) 

A *k + 1 = $ c + Aw k " S 31 Aa k + (S 31 F w K " rK) Av k (169) 


Equations 168 and 169 are an important contribution of the report. In steady state 
the closed-loop aircraft response to deterministic and stochastic disturbances 
can be decomposed into two components. One component, As k , s ^ ows h° w eac ^ 
state in the aircraft is affected by the disturbance states that satisfy a dif- 
ference equation. The other component, Ax^, shows how the white Gaussian noise 
sources are corrupting the system response. Equation 168 leads to a natural frame- 
work for finding the feedback gain K as discussed in the next section. 

* * 

If the solution for Ax^ is substituted into the control system, Au^ can 

be rewritten as 

4 - <- K c * Si - K c „> - 4 (170> 

A \ ' " K C x ^k " K Av k 071) 

In steady-state, the control response decomposes into a component dependent on 
As^ and a zero-mean component, Au^, driven by measurement noise and the zero- 
mean plant response, Ax^. 


(166) 
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D. Optimal Disturbance Suppression: 


The objective of this section is to construct a cost function representing 

desirable objectives for the plant and control system shown in Eqs. 152 to 155, 

then minimize the cost function by finding the feedback gain K. The derivation 

is cumbersome if T and H are non-zero matrices. The matrices r and H are 
w w w w 

assumed to be zero for the derivation in this section. After the optimal gain 

is derived with T and H zero, the optimal gain solution with r and H non-zero 
w w r ° w w 

is presented for completeness. 

If T and H are non-zero, the feedback gain, K, can affect values in the 

matrix and the response of Ax^. The feedback gain cannot affect the response 

of As. . Hence, a useful cost function with F and H zero for disturbance min- 
k w w 

imization is 


J 

e 


|[aS T Au T J 


M 


T 

M R 


Ax 

Au 


+ i trace ](H x S 31 + Hj 1 Q 2 (H x S 31 + HJ } (172) 


( 


If Ax^ were weighted in a quadratic cost function, the As^, component in Ax^ 
would add to the cost. If As^ is neutrally stable or unstable, then quadratically 
weighting Ax^ in a cost function could yield an undefined cost in steady-state. 

Weighting Ax^ and Au^ as shown in Eq.172 yields a finite cost if the closed- 
loop plant matrix is stable, since both signals are zero-mean and do not depend 
on As^. The second part of the cost function in Eq.172 attempts to minimize 
the disturbance response in the system output as defined in Eq. 166 . No restric- 
tion is placed on the first dimension of H and H in Eq. 172 . The cost function 

X s 

represents a tradeoff between optimal stochastic performance and deterministic 
disturbance suppression. 

Substituting Eq. 171 into the cost function produces 
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J = trace (Q. + C 1 K T R K C - M K C - C T K T M T ) X 
e 1^1 x x xx 


+ trace } K 


where 


R K vj 


+ i trace (H x S 31 + H s ) + H g ) 


(173) 


E \ \ " V ’ 


M =x 


X is the covariance of the system error response and is computed using Eq. 169. 
The cost function, J e> can be minimized subject to the equality constraints for 
X and using Lagrange multipliers. The complete cost function is given by 


J = J + J T 
e L 


(175) 


where 


+ trace < (4 S„. + S 01 (-$)+. D -T K C ) Y 

I c 31 31 w w w 


J T = trace J (-X + 4 X 4 T +'W 4-' rKVK T r T + S„. W si, 4 S_. W -W T S^)P T f (176) 
L j c c s 31 31 sx sx 31 \ 


(177) 


4 = 4 - T K C 
c x 


The P and Y matrices are the Lagrange multipliers. 

Necessary conditions for J to have a minimum are 3J/3X = 0, = 0, 

3J/3P = 0, 3J/3Y = 0, and 3J/3K = 0. The five necessary conditions for optimality 
can be expressed as five matrix equations. 


Covariance Equation for x^. 


X = $ X $ + 
c c 


b - s 3i] 


w w 


w w -s:. 

sx s 31 


T T 

i + r k v k r 


(178) 


Star Trajectory Equation 


4 S,. + S,. (-4 ) = T K C - D 

c 31 31 w w w 
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Cost Equation 


T 1 T T 1 

— 

- 


r~ — 

P = <f> P $ + I -C K 

C C 1 X J 

*1 

M 


I 


m t 

R 


-K C 

X 


Complementary Star Trajectory Equation 

Y + Y (-$3) = 2.0 P I W - S_ W , A - H T Q H + H S_. 

C w I sx 31 sj x 2 J_ s x 31 J 

Gain Equation 


( 180 ) 


(181) 


(R + r T P D K (C X C T + V) = (r T P 4 + M T ) X C T + % r T Y (C T + C T ) (182) 

XX X w31x v/ 

A numerical algorithm for finding a gain, K, which satisfies the necessary con- 
ditions is presented next by extending the convergent ouput feedback algorithm 
derived in Ref. 8: 

Choose an initial K so that $ is stable. 

o c 

Choose a scalar a so that J(K^) is decreasing during the following itera- 
tions: 


1. Solve Eq. 179 for using K^. 

2. Solve Eq. 178 for X using and S^ from 1. 

3. Solve Eq.180 for P using K^. 

4. Solve Eq.181 for Y using K^, S^ from 1 and P from 3. 

5. Solve Eq.182 for K using S^ from 1, X from 2, P from 3 and 

Y from 4. 

6. Update the feedback gain 

K i+1 = K i + 01 (K - V (183) 

7. If | J(K^ + ^) - J(K^)| > £ and/or j 1 9j(K_^)/3K^| | > £, where e is 
some convergence criterion, then set i = i+1 and repeat steps 1 
through 6. 

Equations 178 to 181 are Lyapunov equations that can be solved using the 
Bartels-Stewart algorithm. Ref. 19. 
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If H and/or T are non-zero, it is no longer clear that the cost func- 
w w 

tion shown in Eq. 172 represents design objectives. Alternative cost func- 
tions with H and T non-zero is an area of future research. If As, remains 
w w k 

uncontrollable or controlling As^ from Au^ is not desired, as is the case for 
wind shear disturbances, then the cost function shown in Eq. 172 remains attrac- 
tive. The following equations show how the problem is modified if and/or 

T are non-zero: 
w 

First the following matrices are defined': 


r = 

s 

r - 

S 31 

r 

w 

(184) 

$ = 

$ - 

r k t 

■« 

(185) 

c 



X 

r = 

r k 

c - 

H 

(186) 

c 

w 

X 

w 

'v 

$ = 

$ - 

r k 

C 

(187) 

d 

w 

w 

w 

$ a 
cd 

= -9 + S_. 

c 31 

r 

c 

(188) 

dc 

= -*d 

+ r 

c 

S 31 

(189) 


then the five necessary conditions can be shown to reduce to the following 
matrix equations 

Covariance Equation for X 



r EVKT' 
s s 


(190) 


Star Trajectory Generalized Riccati Equation 


$ s„. + s„. (-$,) = r k c - d - s„.rs,; 

c 31 31 d w w 31 c 31 


Cost Equation 


p ■ 4 p * cd ♦[!-<£ x T ] 

*1 

M 


— — 
I 


h t 

R 


-KC ; 
— x 


(191) 


(192) 
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Complementary Star Trajectory Equation 


T T 

$ , Y + Y 
cd dc 


2.0 P W - S„ W | - H T Q, H + H S,. 

I sx 31 s x V I s z 311 

2.0 p $ , x r T - 2.o p r k v k t r T 

cd c s w 


(193) 

(193) 


Gain Equation 

(r + r T p r ) k (c x c T + v) = (r T p ($ - s 01 h ) + m t ) x c t 

ss xx s 31 w x 

+ % r T Y (C T + S* C T ) (194) 

s w 31 x 

The algorithm shown for Eqs. 178 to 182 can be used to solve Eqs . 190 to 194. The 

primary differences are that K must be chosen so that is stable and is 

the solution to a generalized Riccati equation rather than a Lyapunov equation. 

K stabilizes the stochastic term. Ax, when H and T are non-zero, but the total 

w w 

closed-loop system for Ax and As can be unstable or neutrally stable because of 
the way J is chosen. 
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6. DESIGN AND SIMULATION OF CONTROL SYSTEMS 


This chapter presents digital flight control system design parameters 
and nonlinear 6-DOF simulations for four digital control designs. The objective 
of the control system is to cause a small commercial jet (Boeing 737) to capture 
and track a 3 deg glideslope near the terminal area. Enroute to touchdown, the 
SRI wind package. Ref. 31, is used to model a windshear disturbing the ..aircraft 
dynamics . 

A. Design Models 

The four control designs are: 

SR-PIFCGT A. single rate PIFCGT design using a linear aircraft model driven 
by white noise sources. 

MR-PIFCGT A multi-rate PIFCGT design using a linear aircraft model driven 
by white noise sources. 

WS-PIFCGT A single rate PIFCGT design using a linear aircraft model com- 
bined with the windshear model. The feedback gains are computed 
using the algorithm in section 4D with T w and H w nonzero. 

EP-PIFCGT A single rate PIFCGT design using a linear aircraft model combined 
with the windshear model and the energy probe measurement. The 
feedback gains are computed using the algorithm in Section 4D with 

F and H nonzero, 
w w 

The states in the aircraft design model and the output vector used for 
feedback are presented next: 

Integrators 

A W - A5 hk + At [ (A \- A Vk> + 5 - ^m.kU 


A W+1 ’ &5 CAS,k + At (ACAS k - ACAS m,k > 
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Prefilters 


AS = - 10.0 (AS - Aa ) 
z z z 


Aq = - 10.0 (Aq - Aq) 
Complementary Filter 


A 


“ 





— — 

r ~ 1 

— — 

Ah 

= 

0 

1 


Ah 

+ 

0.8 

Ah - Ah + 

0 

■ • 

. A 

Ah 


0 

0 


A 

Ah 


0.32 


1 


SR-PIFCGT At = 0.1 sec 
States 

Ax T = f" Au Aw Aq A9 Ah Aa Aq Ah Ah AEPR AS A£, A£ _ Au Au 

I Z 6 n LAj L 6 


Controls 
Au T = [ A v e AvJ 

Outputs 

AY'*" = |Aa z Aq ACAS A0 Ah Ah A£^ A^^g ^ U t ^ U e] 

MR-PIFCGT At f = 0.05 sec; At =0.1 sec 

I s 

States 

Ax^ = j^Au Aw Aq A0 Ah Aa^ Aq Ah Ah AEPR A6 g AE^ AE^^, Au 
Controls 

Au T - [Av e Av t ] 

Fast Outputs 

■ K A ’ A5 h A5 CAS Au t Au e] 



Slow outputs 

ACAS A0 Ah Ah ] 
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WS-PIFCGT 


At = 0.1 sec 


States 

Ax T = [au Aw Aq A0 Ah Aa Aq Ah Ah AEPR AS A?, A£ Au Au 1 
1 z 0 n LAb l 0 1 

Disturbances 

asT - Kb Aw wJ 

Controls 
Au T - [Av e Av ( ] 

Outputs 

Ay T = |^Aa z Aq ACAS A9 Ah Ah A^ h A?^ Au t Au g J 

EP-PIFCGT At = 0.1 sec 
States 

Ax T = |^AEP 1 AEP 2 Au Aw Aq A0 Ah Aa z Aq Ah Ah AEPR AS £ A? h A? CAg Au t Au g 

Disturbances 
As 1 - [AU„ b AwJ 

Controls 
Au T = [Av e AvJ 

Outputs 

Ay T = [aeP Aa z Aq ACAS A0 Ah Ah A£ h A£ CAg Au t Au £ J 

The aircraft model is trimmed at 69 m/s (135 kt) descending straight flight 
using a glideslope of 3 deg. Flaps are set to 30 deg. The weight is 85000 lbs. 
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The five states for the longitudinal aircraft dynamics in body axis are 

the body x-axis velocity, Au, the body z-axis velocity, Aw, pitch rate, Aq, 

pitch Euler angle, A0 and vertical height Az. The prefilter states, Aa z and Aq, 

model a continuous-time prefilter used to suppress aliasing errors in the body 

mounted z-axis accelerometer measurement for a and the rate gyro measurement 

z 

for q. The noisy height measurement, h, which can come from a barometric altimeter 

or the Microwave Landing System (MLS), is complementary filtered with the gravity 

* 

A /\ 

corrected output of the accelerometer to yield smooth estimates for h and h. The 

two states, AEPR and A6 g , model the actuators dynamics for throttle and elevator 

respectively, and are discussed in Ref. 32. The type 1 property for the control 

system is obtained using the integrator states L and £ . The integrator 

n laj 

state, A'£ , integrates the height command error plus 5 times the height rate 
error. The integrator state, A£ , integrates the CAS command error where 
CAS is a calibrated airspeed measurement. The states Au^ and Au^ are the per- 
turbation aircraft control states placed in the state vector because of control 
rate weighting. 

The multi-rate control design breaks the measurements into a fast-sample- 
rate group and a slow-sample-rate group. The fast rate is chosen to be 20 
samples/sec while the slow rate is 10 samples/sec. The integrator states and 
control position states must be sampled at the fast rate. The two inner-loop 
measurements, Aa^ and Aq, are also sampled at the fast rate. The slow measure- 
ments are airspeed, pitch angle and the output of the complementary filters. The 
complementary filter model in $ is not modeled at the slow rate, however. Model- 
ing a slow rate complementary filter in $ would cause $ to become a periodic 
matrix. The complementary filter is propagated and updated at 20 iterations/sec 
in the nonlinear simulation. The WS-PIFCGT design uses the wind shear disturbance 

dynamics in the design model. The nominal values used for u and w are -0.1 m/ 

z z 

sec/m and 0.02 m/sec/m, respectively. 
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The EP-PIFCGT design includes the energy probe measurement in the WS-PIFCGT 
measurement vector. Two states, AEP^ and AEP^, are added to the state vector to 
model the sensor dynamics described in Refs. 10 and 11. 

The longitudinal command model propagated in the simulation uses a double 
integration of vertical acceleration to generate the vertical path trajectory. 
The CAS trajectory is constant, 



(195) 


The commands used for "h in the nonlinear simulation are 

m 


If 

h 

> -8.3 

and 

h 

> 200 

Then 

h = -1.5 Else 

(196) 


m 



m 



m 


If 

h 

< 150 

and 

h 

< 0 

Then 

• • *t 

h =1.5 Else h =0.0 

(197) 


m m m m 


The command model descends at a constant rate of descent (8.3 ft/sec) until 
150 ft is reached. The descent is changed to straight and level flight after h^ 
reaches 150 ft. The change to straight and level flight eventually puts the 
command model height state at 122ft. . The aircraft and command model are initial- 
ized to 700ft straight and level flight at the beginning of the simulation. The 
command model used to design the command generator tracker feedforward gains only 

t 

uses h , h and CAS as discussed in Ref. 3. 
m m m 

The design of the control system is achieved by choosing the diagonal ele- 
ments in the matrices Q, R, X, W and V iteratively until the control system has 
desired closed-loop properties. Adjustments are made to the elements until one 
set of numbers, shown in Table 1, can be used in all four design conditions . 
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B. Control Design Properties 


Six evaluation procedures are used to determine the feedback and tracking 
properties of the control designs. The six evaluations are the feedback gain ele- 
ments, singular value plots, Ref. 33, eigenvalue plots, Ref. 33, Bode plots for 
each control loop, closed-loop eigenvalues and nonlinear simulations. 

The feedback gain and feedforward gain matrices are shown in Table 2. Com- 
paring the single rate and multi-rate designs, the gain elements almost all in- 
creased in value an average of 25% in the elevator loop. The elevator loop band- 
width increased in the multi-rate design as shown by the -9.37 control gain. The 
multi-rate design had less of an effect on the throttle loop which has a much 
lower bandwidth. The two wind shear design gain elements change significantly 
from the case when no wind shear is modeled in three of the measurement loops: 

A A 


a z , h and One of the closed-loop poles is significantly altered (de- 

stabilized) by the wind shear dynamics and the increased gain values on these 
three measurements appear to improve system robustness to wind shear distur- 
bances. The throttle loop bandwidth increases in the presence of the wind shear 
dynamics as evidence by the -3.63 and -4.9 5 T control gains. 

The closed-loop eigenvalues for the designs are obtained from the discrete 
system linear matrices. Each eigenvalue inside the unit circle in the z-domain 
is mapped to the left-half plane using 

je 


X z = a + j b = me J 


(198) 


1-1-, } _ m , . 9 

s At n z At J At 


(199) 


The three single rate designs are evaluated with and without the wind shear 
dynamics. The SR-PIFCGT closed-loop eigenvalues shown in Table 3 . have adequate 
damping and stability without the wind shear system. The wind shear dynamics 
introduce two poles into the system; one is neutrally stable while the other is 
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slightly unstable. Windshear as a phenomenon only occurs for a specific period 
of time. The control system's objective is to maintain control over the aircraft 
dynamics and track the command model. The primary change in the closed-loop 
eigenvalues is that the pole most associated with the phugoid motion (oj^ = 0.35, 

£ = 0.49) changes damping to £ = 0.19, i.e., the phugoid motion destablizes. 
Increasing the nominal values for the wind shear gradients further destabilizes 
the phugoid mode. 

The new capability provided by the theoretical derivation in Section 5D allows 
the design plant model to contain the unstable coupled windshear dynamics. The 
WS-PIFCGT closed-loop eigenvalues are presented in Table 4.' The closed-loop 
model without wind shear has good damping ratios on complex modes (£ > 0.5) with 
the windshear model, the phugoid mode destabilizes but has better damping than 
the mode in Table 3 (5 = 0.33 verses £ = 0.19). The- energy probe further aids 
in stabilizing the phugoid mode as shown in Table 5 (£ = 0.40). 

The closed-loop mapped eigenvalues for the multi-rate design are shown in 
Table 6. Despite a significant change in feedback gains between single and multi- 
rate designs, the closed-loop eigenvalues for the single rate and multi-rate 
designs are similiar. The multi-rate and single rate designs use the same qua- 
dratic weights. The extra three poles in Table 6 are caused by the hold cir- 
cuits in the system model. 

The frequency domain methods for evaluating the single rate control designs 
result in the plots shown in Figs. 1 to 6 and Table 7. Frequency domain an- 
alysis of the multi-rate design is not within the scope of this effort. The 
minimum singular value of the return difference matrix, as discussed in Ref. 33, 
is a conservative indication of closed-loop plant robustness to unstructured 
perturbations in the plant dynamics. A small singular value means there is a 
small plant perturbation (physically unrealizable perturbations are allowed) 
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that can destabilize the plant. The three singular value plots all indicate 
potentially poor designs. Reference 33 investigated singular values and demon- 
strated that eigenvalue analysis of the return difference matrix, although a 
not guranteed indication of robustness, provides a clear picture of the poss- 
ible robustness boundary. The eigenvalues of the return difference matrix are 
bounded below by the minimum singular value and above by the maximum singular 
value. The true multivariable plant robustness should be bounded between the 
minimum singular value and the minimum eigenvalue of the return difference matrix. 

The eigenvalue plots in Figs. 1 to 3 indicate a much more robust design then 
the singular value plots and compare favorably to the information provided by 
the Bode plots. The smallest eigenvalue is near to = 3.5 rad/sec in all three 
design. Specific values are shown in Table 7. The frequency for the smallest 
eigenvalue is near the frequency for the short period modes shown in Tables 3 to 5. 
Another effect evident from the eigenvalue and singular value plots in Fig. 3 is 
that the energy probe sensor increases plant robustness in the low frequency 
region, u)<l, where plant alterations due to windshear are likely to occur... 

The Bode plots in Figs. 4 to 6 provide useful information concerning the 
bandwidth and robustness of each loop, individually. The bandwidth (the frequency 
at which the gain remains below -6 db) for throttle is between 0.5 rads/sec and 
0.25 rads/sec. The designs using the windshear dynamics have lower throttle 
bandwidth. The bandwidth for the elevator is between 2 and 3 rads/sec. All phase 
margins are better than 65 degs and gain margins are better than -10 db. 

C. Nonlinear Simlations 

Nonlinear six degree-of-freedom simulations of the control systems for a 
glideslope capture and track are shown in Figs. 7 to 9 . Windshear, gust dis- 
turbances and measurement noise are inactive in the simulation. The captures are 
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smooth with little overshoot both at the beginning and end of the simulation where 
the aircraft returns to straight and level flight. The best performance is pro- 
vided by the multi-rate design which has smoother surface motion and better CAS 
control . 

The next four simulations in Figs. 10 to 13 show the four control designs 
capturing and tracking the glideslope in the presence of windshear. The longi- 
tudinal, horizontal and vertical shears are shown in Fig. 16 and correspond to 
the windshear in the SRI package known as the Philadelphia/Allegheny windshear. 

The. aircraft encounters an increasing headwind. The headwind decreases grad- 
ually then dramatically into a tailwind. The tailwind changes back to zero. Coin- 
cident with the change from a headwind to a tailwind, the aircraft encounters a 
downward vertical wind of 20kt. The large drop in airspeed causes the control 
systems to briefly saturate throttle (at 60 deg) and pitch the aircraft up to 
8 deg. 

The best windshear response is provided by the WS-PIFCGT design. The air- 
speed deviations are smaller and the aircraft pitches up the least (6 deg) . The 
WS-PIFCGT design also maintains the largest separation between the aircraft and 
the ground in the critical time period between 75 and 85 secs into the simula- 
tion. The MR-PIFCGT design has the largest deviations in airspeed and height 
during the windshear. The EP-PIFCGT design using the energy probe sensor for 
feedback encountered small elevator and throttle oscillations when throttle 
surpassed 30 degrees. 

A 1 kt gust is added to the nonlinear simulations in Figs. 14 and 15. The 
effect of the gusts is primarily evident in the throttle response. The tracking 
is relatively unaffected by the gust response. 
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7 . SUMMARY AND RECOMMENDATIONS 


A summary of the results in this report are: 

• Optimal output feedback is successfully used to design a single rate 
control system using the proportional integral filter structure. 

• The optimal output feedback approach and optimal output feedback 
algorithm are extended to use the hierarchical control design approach 
with minimal change to the output feedback design model matrices. Un- 
wanted states and observations encountered as each loop is designed are 
made uncontrollable and unobservable and stabilized if the mode is not 
stable. 

• An optimal multi-rate output feedback synthesis problem is formulated and 
necessary conditions for an optimal solution are derived. An algorithm 
to solve for the constant gain solution is presented. The algorithm is 
programmed and designs are made comparing the multi-rate and single rate 
designs. The PIF structure is extended to use multi-rate output feed- 
back. The PIF feedforward gain is shown to be periodic in time. 

• The command generator tracker is extended to accommodate disturbances 
which are affected by the plant dynamics. The extension causes one of 
the Lyapunov equations encountered in the usual command generator tracker 
to become a generalized Riccati equation. A globally convergent algorithm 
for solving the generalized Riccati equation is developed. 

• A new and novel approach is developed to accommodate windshear in the 
presence of noise. The command generator tracker is extended to the 
stochastic environment. The steady-state plant response is decomposed 
into the stochastic and deterministic components. A novel cost function 
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is constructed which suppresses the effect of disturbances, stabilizes the 
stochastic plant response and does not require a measurement of the dis- 
turbance state. The disturbance and plant dynamics can be coupled. An 
algorithm which computes the (local) optimal gain is developed, has been 
programmed and is used to determine a PIFCGT control system. 

• Design and simulations showed that the multi-rate design (which uses 
faster measurements on some sensors) performs best for the quadratic 
weights chosen when no disturbances are present. The design using the 
new disturbance accommodation synthesis approach performs best when the 
aircraft is simulated flying through a windshear. The windshear design 
using the energy probe also performs well except for small oscillations 
that occur during the severe part of the windshear simulations. 

Recommendations 

m The command generator tracker can be extended to the multi-rate case. 

• The multi-rate command generator tracker would allow the computation of 
the multi-rate frequency response by extending the approach developed 
in Ref. 30. 

•’ The multirate problem can be extended to include the case of many 
models with one controller as discussed in Ref. 17. 

• The solution in Section 5B has application beyond disturbance suppress- 
ion. The solution in Section 5B solves a feedback decoupling problem. 
Consider the lateral dynamics of the aircraft where it is desirable to 
decouple the roll and sideslip response. Partioning the plant states 
into two groups and identifing one group as Ax and the other group as 
As, the feedforward solution yields a control system for decoupling the 
HAx response from the As states. HAx can be sideslip and As can be the 
roll state. Any rolling motion would not affect sideslip. 
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APPENDIX A 


A Convergent Algorithm for a Generalized Riccati Equation : 

A generalized Riccati equation has the following form: 

X = AXB + AXCX + D (200) 

where X is an unknown matrix. If A were invertible then Eq. 200 could be 
rewritten as 


0=FX+XB+XCX+G (201) 

which resembles the continuous- time Riccati. equation obtained in the optimal 
regulator problem. 

Define the function f(X) as 

f(X) = AXB + AXCX + D (202) 

and define the direction, d, at X as 

d=AdB+AdC+AXCd+ (f(X)-X) (203) 

Define the matrix, E, as 

E = (I - A X C) -1 (204) 

the direction, d, satisfies the Lyapunov equation 

d = E A d (B + C)+E (f(X)-X) (205) 


The globally convergent algorithm for computing X is 

1. Choose a starting value for X, usually X^ = 0. 

2. Compute d in Eq. 205 using X . 

K- K. 

3. Choose according to 



min 


1 . 0 , ^ 


II ^V-Xkll j 

II A d k C d k || ) 


C206) 
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(207) 


4. Update X^ using 


X k+1 =X k +a k d k 


5. If ||;f (X k ) - X ) || is small, stop, otherwise increment k and go back 


to 1. 
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TABLE 1. DESIGN PARAMETERS FOR THE CONTROL DESIGNS 
(UNITS ARE FT. AND DEG) 
















CONTROL 


STATE FEEDBACK GAINS 


CONTROL DESIGN 


E 

P 


A 



q 


CAS 0 h 


h 


V T 

- 

- 6.33 

- 0.566 

- 2.18 

- 0.47 

- 0.57 

- 2.69 

1 

v e 


1.65 

4.23 

1.91 

6.69 

0.76 

1.72 

v„ 


- 9.38 

- 0.37 

- 2.10 

0.171 

- 0.923 

- 3.10 

1 

v e 


5.14 

4.80 

1.91 

6.88 

1.06 

4.09 

V T 

- 1.35 

- 9.86 

- 0.90 

- 2.05 

1.14 

- 1.06 

- 4.97 

i 

V 

e 

1.14 

4.69 

4.08 

2.67 

7.56 

1.08 

4.92 

v T 


- 6.92 

- 0.41 

- 2.24 

- 0.367 

- 0.628 

- 2.86 

l 

V 

- 

2.56 

5.5 

2.39 

8.22 

0.99 

2.27 


SINGLE RATE 
DESIGN 


SINGLE RATE 
WINDSHEAR DESIGN 


SINGLE RATE 
WINDSHEAR DESIGN WITH 
ENERGY PROBE 


MULTI-RATE 

DESIGN 








TABLE 2. FEEDBACK GAINS FOR THE CONTROL SYSTEMS (CONTINUED) 


CONTROL 

INTEGRATOR GAINS 

CONTROL DESIGN 


CONTROL 

CONTROL 

GAINS 

CONTROL DESIGN 





6 t 

6 

e 


V T 

-0.075 

-0.11 

SINGLE RATE 


V T 

-2.59 

-1.0 

SINGLE RATE 

V 

e 

0.064 

0.003 

DESIGN 


V 

e 

1.31 

1.76 

DESIGN 

v 

-0.086 

-0.184 

SINGLE RATE 


V T 

-3.63 

-0.75 

SINGLE RATE 

V 

0.099 

0.060 

WINDSHEAR DESIGN 


V 

2.37 

-8.53 

WINDSHEAR DESIGN 

e 





e 




V T 

-0.196 

-0.26 

SINGLE RATE 


v 

-4.9 

0.174 

SINGLE RATE 

V_ 

UNAVAILABLE 

WINDSHEAR DESIGN WITH 


v e 

3.75 

-7.1 

WINDSHEAR DESIGN WITH 




ENERGY PROBE 





ENERGY PROBE 

V T 

-0.077 

-0.11 

MULTI-RATE 


V T 

-2.64 

-0.631 

MULTI-RATE 

V 

0.083 

0.0078 

DESIGN 


V 

1.67 

-9.37 

DESIGN 

e 





e 





w 



















TABLE 2. FEEDBACK GAINS FOR THE CONTROL SYSTEMS (CONCLUDED) 


CONTROL 

FEEDFORWARD 

GAIN 

CONTROL DESIGN 


h 

m 

V 

m 


V T 

' 5.4 

1.4 

SINGLE RATE 

V 

e 

-3.4 

-0.039 

DESIGN 

V T 

7.73 

0.785 

SINGLE RATE 

V 

e 

-6.59 

0.157 

WINDSHEAR DESIGN 

V T 

7.89 

0.373 

SINGLE RATE 

V 

e 

-6.59 

0.157 

WINDSHEAR DESIGN WITH 
ENERGY PROBE 

MULTI-RATE DESIGN 

v 

-0.064 

-0.016 

-K^B 

v e 

0.032 

-0.0033 

V T 

2.95 

2.19 

-K CA n 

i 

V 

e 

-4.21 

-1.34 

s s 12 

v T 

1.65 

-1.26 

-K^A -K-C..A, „-K A„„ 

T 

V 

e 

0.75 

1.26 

£ f f 12 u 22 
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TABLE 3. CLOSED LOOP MAPPED EIGENVALUES FOR THE 
SINGLE RATE DESIGN (SR-PIFCGT) 


DESIGN CONDITION 
WITHOUT WINDSHEAR 

DESIGN CONDITION 
WITH WINDSHEAR 

IDENTIFIABLE MODE 

REAL 

a) 

_.n 

(rad/sec) 

G 

REAL 

w 

n 

(rad/sec) 

G 



21.3 

0.84 


21.3 

0.84 


-10.0 



-10.0 




- 4.7 

2.60 

0.58 

- 4.7 

2.60 

0.58 

SHORT PERIOD 


0.84 

0.98 


0.82 

0.99 



0.52 

0.68 


0.63 

0.69 



0.35 

0.49 


0.369 

0.19 

PHUGOID 


0.079 

0.64 


Q. 063 

0.90 


-0.08 



-0.13 







0.0 







0.04 







TABLE 4. CLOSED- LOOP MAPPED EIGENVALUES FOR THE 

SINGLE RATE DESIGN WITH WINDSHEAR (WS-PIFCGT) 


IDENTIFIABLE MODE 








TABLE 5. CLOSED-LOOP MAPPED EIGENVALUES FOR THE SINGLE RATE 
WINDSHEAR DESIGN WITH THE ENERGY PROBE (EP-PIFCGT) 


DESIGN CONDITION 
WITHOUT WINDSHEAR 

DESIGN CONDITION 
WITHOUT WINDSHEAR 

IDENTIFIABLE MODE 

REAL 03 Z 

REAL 03 4 


n 

n 


(rad/ sec) 

(rad/ sec) 


22.13 0.71 

22.13 0.71 


-10.0 

-10.0 


- 4.5 

- 4.5 


3.74 0.56 

3.74 0.57 

SHORT PERIOD 

1.65 0.89 

1.64 0.91 


0.90 0.55 

0.88 0.55 



0.57 0.77 


0.21 0.74 

0.37 0.40 

PHUGOID 

0.08 0.41 




0.072 0.57 


-0.13 

-0.08 



0.0 



0.039 
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TABLE 6. CLOSED-LOOP MAPPED EIGENVALUES FOR THE MULTI-RATE 
CONTROL DESIGN 


DESIGN 

WITHOUT 

CONDITION 
WIND SHEAR 


REAL 

0) 

n 

(rad/sec) 

C 


370 

0.99 


366 

0.99 


325 

0.99 

-10.0 
- 5.2 

20.6 

0.96 


2.6 

0.57 


0.82 

0.95 


0.52 

0.69 


0.34 

0.50 

-0.075 

0.080 

0.66 


IDENTIFIABLE MODE 


SHORT PERIOD 


PHUGOID 








TABLE 7. FREQUENCY DOMAIN PROPERTIES OF THE CONTROL DESIGNS 

(to is in rad/sec) 


CONTROL 

DESIGN 

MINIMUM 
SINGULAR VALUE 

MINIMUM 

EIGENVALUE 

THROTTLE 

ELEVATOR 

GAIN 

MARGIN 

PHASE 

MARGIN 

GAIN 

MARGIN 

PHASE 

MARGIN 


db 

10 

db 

10 

db 

<0 

deg 

to 

db 

to 

deg 

to 

SR-PIFCGT 

-15.07 

0.4 

-2.8 

3.3 

-18.0 

1.4 

84 

0.15 

-13.0 

4.6 

80 

0.45 

WS-PIFCGT 

-12.8 

0.9 

-3.7 

3.3 

-20.0 

1.7 

70 

0.1 

-40.0 

4.3 

80 

0.63 

EP-PIFCGT 

-14.1 

0.1 

-3.5 

3.9 

-20.8 

1.7 

84 

0.15 

-10.0 

4.7 

69 

0.57 






























SINGULAR VALUES (dB) EIGENVALUES (dB) 


a) Eigenvalues (dB) , | |-High value, Q-Low value 
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FIGURE 2. SINGULAR VALUE AND EIGENVALUE PLOTS FOR THE SINGLE RATE WINDSHEAR 
DESIGN (WS-PIFCGT). 







SINGULAR VALUES (dB) EIGENVALUES (dB) 


a) Eigenvalues (dB), £”]-High value, Q-Low value 


FIGURE 3. 



b) Singular values (dB) , Q-High value, | |-Low value 



SINGULAR VALUE AND EIGENVALUE PLOTS FOR THE SINGLE RATE WINDSHEAR 
DESIGN WITH THE ENERGY PROBE (EP-PIFCGT) 
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GAIN (dB) GAIN (dB) 
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FIGURE 6. 

BODE PLOTS FOR 
ENERGY PROBE. 

THE SINGLE RATE WINDSHEAR DESIGN 
(EP-PIFCGT) 

WITH THE 


PHASE (deg) PHASE (deg) 











LONG CONTROLS (DEG) EPR THR0TLE/20 DEGREES 





FIGURE 7. SR-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH NO WINDSHEAR AND NO GUSTS 
(CONCLUDED) 
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FIGURE 8. MR-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH NO WINDSHEAR AND NO GUSTS 
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LONG CONTROLS ( DEG) 


0 10 20 30 >10 SO 60 70 80 30 100 

TIME ( SEC) 



FIGURE 8. MR-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH NO WINDSHEAR AND NO GUSTS 
(CONCLUDED) 
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FIGURE 9. WS-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH NO WINDSHEAR AND NO GUSTS 
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LONG CONTROLS (DEG) EPR THR0TLE/20 DEGREES 





FIGURE 9. WS-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH NO WINDSHEAR AND NO GUSTS 
(CONCLUDED) 
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CAS (M/S) , HEIGHT (M) 





FIGURE 10. SR-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WENDSHEAR AND NO GUSTS 
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LONG CONTROLS ( DEG) EPR THR0TLE/20 





FIGURE 10. SR-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WINDSHEAR AND NO GUSTS 
(CONCLUDED) 
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FIGURE 11. MR-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WINDSHEAR AND NO GUSTS 
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FIGURE 11. MR-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WINDSHEAR AND NO GUSTS 
(CONCLUDED) 
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HDOT (M/S) 


TIME (SEC) 



TIME ( SEC) 


FIGURE 12. WS-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WINDSHEAR AND NO GUSTS 



LONG CONTROLS ( DEG) EPR THR0TLE/20 





FIGURE 12. WS-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WINDSHEAR AND NO GUSTS 
(CONCLUDED) 
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LONG CONTROLS ( DEG) EPR THR0TLE/20 DEGREES 





FIGURE 13. EP-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WINDSHEAR AND NO GUSTS 
(CONCLUDED) 
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FIGURE 14. SR-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WINDSHEAR AND A lkt GUST 
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LONG CONTROLS (DEG) EPR THR0TLE/20 DEGREES 





FIGURE 14. SR-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WINDSHEAR AND A lkt GUSTS 
(CONCLUDED) 
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125 





FIGURE 15. MR PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WINDSHEAR AND A lkt GUST 
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LONG 



FIGURE 15. MR-PIFCGT GLIDESLOPE CAPTURE AND TRACK WITH WINDSHEAR AND A lkt GUST 
(CONCLUDED) 
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